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A generic mechanism of collapse in the Gross-Pitaevskii 
equation with attractive interparticle interactions is gained 
by reformulating this equation as Newton's equation of mo- 
tion for a system of particles with a constraint. 'Quantum 
pressure' effects give rise to formation of a potential barrier 
around the emerging singularity, which prevents a fraction 
of the particles from falling into the singularity. For reason- 
able initial widths of the condensate, the fraction of collapsing 
particles, which are thereby removed from the condensate, is 
found to be a 'universal' number ~ 0.7. 



The Bosc-Einstein condensate formed in a magneto- 
optically trapped vapor of cold alkali atoms, can be quite 
accurately described by the Gross-Pitaevskii (GP) equa- 
tion@ 



2m 



A* -ff*|*| 2 - V(x)V = 0, 



(1) 



where ^(x, t) is the wave function of the condensate, 
the external potential V(x) models the wall-less confine- 
ment (the trap), m is the mass of an individual atom, 
g — ATrh 2 m~ 1 a s is the effective interaction strength (the 
'coupling constant') with a s the scattering length, and 

A = 22 i ~§^z i s the Laplace operator. A convenient 
as well as practical choice for the confining trap is the 
paraboloidal potential V — y Yli ljJ 1 x %' ■ 

It has been known Jl]-|7j for some while already that 
the Bose-Einstein condensate of atoms with a negative 
scattering length becomes unstable when the number of 
atoms exceeds some critical value N c . These instabilities 
have recently been analyzed ^ experimentally. To this 
end the most interesting prospects are provided by the 
application || of Feshbach resonance by which one can 
actually tune the effective scattering length with an ex- 
ternal magnetic field, and thereby e.g. sweep the coupling 
constant from positive to negative. This technique was 
very recently used § on a condensed gas of 85 Rb atoms 
such that a collapse and the consecutive 'explosion' of a 
condensate of 85 Rb atoms was created. 

Triggered by these recent experiments, a possible 
mechanism through elastic collisions for the ejection from 
the condensate of atoms with negative scattering length 
was discussed in JlOfl . In this Letter we show however 



that the mean-field theory (i.e. the GP equation) of the 
condensate already allows for a description of 'partial col- 
lapse', and leads to an accurate estimate for the number 
of particles that remain in the condensate after the ex- 
plosion. The fraction of the condensate thrown away in 
the collapse (i.e. the collapsing part of the condensate), 
is always about 70% provided that the width of the con- 
densate before the change of sign of the scattering length 
is close to that of the ground state. This was true e.g. for 
the experiment reported in [^) . We find that the collaps- 
ing fraction is relatively insensitive to the total number of 
particles in the condensate, and that this fraction tends 
to 100% only in the limit of very large width of the con- 
densate. 

Collapse phenomena are very well known in nonlinear 
optics and plasma physics, and they have there a long 
history extending over the last thirty years (see e.g. re- 
view [jll] and references therein). Following this tradi- 
tion, the collapse in the GP system with a negative scat- 
tering length was discussed in a similar fashion in [^2[-pT| . 

Within this tradition, the analysis of singularity for- 
mation is mainly focussed on the self-similar behaviour 
of the solutions in the vicinity of the singularity. We have 
recently shown Jig ], however, that a true self-similar so- 
lution of the GP system is only possible in two space 
dimensions. This means that a self-similar description 
of collapse can only be used with certain reservations. 
In this Letter we propose a completely different descrip- 
tion for the singularity formation, which does not rely 
on self-similarity of the solution. In a sense the starting 
point of our analysis is the concept of 'partial collapse', 
now of the Bose-Einstein condensate. A suggestion about 
a singularity restricted to a localized domain inside the 
condensate, was already made in Here we elabo- 

rate further on this suggestion, and formulate it within a 
completely different framework from that used in . 

Our approach is based on the realization that the GP 
equation Eq. (Q) can be interpreted as Newton's equa- 
tions of motion for (fictitious) particles with a certain 
constraint. The dynamics of GP collapse can then easily 
be followed by following the trajectories of these parti- 
cles, under the influence of an effective potential due to 
'pressure' effect that arise from the spatially varying con- 
densate density. By introducing an auxiliary field r)(x,t), 
and its inverse r(T], t), we can show (see below) that these 
trajectories satisfy a system of ordinary differential equa- 
tions, 



1 



dr.i 1 d(j) 

-sr = —-^—i r (V,t),t); 

at m axi 



n(ri,0) = in, 



(2) 



with i = 1,2,3 for d=3 dimensions. Here is a phase 
defined by the representation <3> = ^~pe l ^/ h of the conden- 
sate field. The fields ri(r],t) can now be interpreted as 
trajectories of some particles with initial points r\i. We 
describe in the following how the GP equation Eq. (|l|) 
can be transformed into the equations of motion for this 
system of particles. 

It is plain that, by using \E" = ^/pe^^, Eq. (|l|) is equiv- 
alent to a coupled set of equations for p and 0, 



pt H div(pV0) = 

m 



(3) 



+ V + gp-^p-^Ap^ = 0. (4) 



2 m 



Let us now introduce an auxiliary field rj(x,t), which 
is assumed to satisfy 



at m 
If we define a Jacobian 

J(x, t) = det J™, J™ = 
it then follows from Eq. (^) that 



(5) 



dx n 



^1 + I V In JV<£ + - = 0. 
at m m 

From this Eq. (||) we can easily deduce that 
p(x,t) = J(x,t)po(rj) 



(6) 



(7) 



is a solution of the continuity equation Eq. (|) with po 
an arbitrary function. We now define the field r(r),t) by 
requiring that the function r(rj,t) is the inverse of the 
function r\{x,t). With this definition Eq. (||) and Eq. (^) 
become equivalent. Notice that the relation Eq. (Q) pro- 
vides the anticipated representation of the condensate 
density in terms of particle trajectories. 

Thus far we have not used Eq. (Q). By differentiating 
Eq. (0) with respect to t, we find that 



dt 2 



(8) 



Using now Eq. (|), this Eq. (|) can be expressed in the 
familiar form of Newton's equations of motion for a set 
of particles, 



'' dt 2 



dV 



dU 

dxi ' 



(9) 



in which the 'pressure potential' U(x,t) is 



U(x,t) 



2m 



-l/2 Aj0 l/2 



gp- 



(10) 



The 'particles' of this system thus feel, in addition to 
the confining trap potential, an effective potential that 
arises from the spatially and temporally varying conden- 
sate density p(x,t). The second term in the right side of 
Eq. (^) can be recognized as the conventional pressure 
term, while the first term there arises from the kinetic 
energy part of the GP equation, and can be considered 
as a quantum pressure term. As will be explained be- 
low, it is the interplay between these two effects in the 
pressure potential which determines the dynamics of the 
collapse. 

As can be readily seen, the result Eq. can also be 
written in the form 



<^ I P{r,t) = po(v)- 



(11) 



This constraint can now be used to exclude the parti- 
cle density p from Eq. (^). This means that the system 
of equations Eqs. (p[]Tl|) is self consistent, and fully de- 
scribes the dynamics of the condensate. We need however 
to specify first the initial conditions for this system, and 
a natural choice for Eq. (0) is 



(12) 



with v — ^V0|t = o- The function po can now be inter- 
preted as the initial density of the particles, and the phase 
(j> can be solved from Eq. (||) up to a time-dependent fac- 
tor. 

Notice that there is a relation between (f> and p. Indeed, 
from Eq. (|[) we find that 

det f^ =e ij:w, 

Hence, we can express p in the form 

p(x, t) = e"» /o WMM* poMx, t)), (13) 

where x(t, t 1 ) = r(r)(x, t), t'). 

It is evident that a possible collapse of the condensate 
is related to the singular points of the solution, i.e. to 
the points where the function r(rj,t) is not invertible and 
the determinant det(J^) vanishes. At these points the 
density becomes infinite. It seems natural to assume that 
if collapse takes place, it happens at a finite time t = t*, 
i.e. p{r Sl t*) = oo. Then, from Eq. (fTl|), we can conclude 
that for any r s there is a set fi s such that 

hm ri(r},t) = r si 

for all r/ G £l s - In other words, all particle trajectories 
that begin at any f|£!l s are converged at one point r s . 
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Suppose that the set of singular points, {r s }, consists 
of one isolated point while f2 s is a domain. We can then 
define the number of particles involved in the collapse as 



(14) 



This number can thus be interpreted as the number of 
particles that are ejected from the condensate when it 
is collapsed. Notice that the domain £l s can also coin- 
cide with the whole space. In this case N e = N, the 
total number of atoms in the condensate. This is true in 
particular for the self-similar solutions, when the whole 
condensate collapses to a ^-function centered in the trap, 
and completely disappears apart from the very particu- 
lar borderline case of an oscillating solution |15|. The 
self-similar solutions can only exist |l5| in two space di- 
mensions. It is important to emphasize that, in general, 
N e < N. It can be shown, e.g. , that if po tends to a 
Gaussian distribution as \rj\ — > oo, then N e < N. This 
means that A^ is an essential characteristics of the col- 
lapse of attractive condensates. 

We can explain the partial collapse qualitatively as fol- 
lows. If, initially, the width of the condensate is greater 
than the characteristic oscillator length ano, both the 
quantum pressure and the conventional pressure induce 
the condensate to decrease its width. However, when the 
width becomes equal or less than clho, the two pressure 
terms in the energy have opposite effects: quantum pres- 
sure favours an expanding condensate (as a result of an 
'uncertainty principle'), while the negative conventional 
pressure term favours a collapsing condensate. Provided 
that the number of particles is greater than a critical 
value, the negative pressure dominates over the quantum 
pressure in a small domain centered at a point r s . Out- 
side this domain it is the quantum pressure that dom- 
inates, and prevents the outside particles from entering 
this domain. 

In order to estimate the value of N e , we use here a 
Gaussian trial wave function |l7|-|20|] . It gives an approx- 
imate solution for the density, which can be expressed in 
the form 



p{x,t) 



N 



in which the functions Tj satisfy |21|- 

2 H 2 1 gN I i 

r. i +CL) 4 T l -- 







(15) 



(16) 



2 r? (2^) 3 / 2 m Uk Tk n 

with Tj(0) = ctj and f (0) = 0. Here <Xj is the initial width 
of the condensate. 

Substituting the solution Eq. (|l5|) in Eq. ([}]), we find 
that 



o 2 



dV r 



eff 



dt 2 

V eff (r,t) 



dr 
2 ^ 



(17) 



-p{r,t). 



Here V e // is an effective time-dependent potential that 
determines the behaviour of particle trajectories in the 
vicinity of the center of the trap. 

For simplicity, we consider here in more detail the 
symmetric trap: LUi — lu and ai = a. Hereafter we 
measure the time in units of 1/uj, and the distance in 
units of ano = y/h/moj, the oscillator length. Now 
Eqs. Q,© take the form 



fr _ 
dt 2 ~ 

Veff = 



eff 



dr 
1 - 



(18) 



and 



dt 2 



(19) 



, and the initial conditions are given 
n,r(Q) = 0,t(0) = a,r(Q) = 0. Notice that if 
.9o > 3ocr, where go cr ~ 0.671 U, the solution of Eq. (|l|) 
has a zero at a finite time t = i*. 
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FIG. 1. The effective potential V e // at different times 
t < t t . The initially almost paraboloidal potential develops a 
growing barrier at the final stages of the collapse. 

Let us assume that the system has a > 1 and begins 
to collapse. Initially, V e // can be well approximated by 
a paraboloidal trap. The particles are drifted towards 
the center and the condensate width decreases. When 
this width is reduced to the oscillator length, a potential 



barrier is however formed at a distance r = 
the center, 



from 



850 



<Jtt(t - T 5 ) ' 

preventing the outside particles from falling into the 
emerging singularity. Figure 1 illustrates the formation 
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of this potential barrier in the limit when the width of 
the condensate r goes to zero. 

We solved the system Eqs. (|I8|) , (|l9|) numerically. The 
collapse time was obtained from the condition^ r(i*) < 
0.1. Equation ( fL8[ ) was solved in the time interval [0, i*], 
while the point r\ was supposed to belong to the domain 
fi s provided that r(i*) < r ma x(t*). 

The conditions of the experiment reported in || were 
modelled such that the initial width of the condensate 
was derived as a stationary solution of Eq.(|l9|) with go 
replaced by — 41.67 go- Notice that in the experiment jjj 
the initial positive scattering length was 2500ao (ao is 
the Bohr radius), and this scattering length was in 0.5ms 
changed to — 60ao. The number —41.67 is thus the ratio 
of these two scattering lengths. We varied the parame- 
ter go in a large range, from 1 to 300. In all cases the 
fraction of the particles ejected from the condensate as 
a result of collapse, was about 70%. We also found that 
the fraction of ejected particles increased above 70% only 
when the initial width of the condensate was increased 
well above that of the ground state. In the limit of very 
large widths the fraction approached 100% of the total 
number of particles. 

In summary, in this Letter we have explained the 
generic mechanism of collapse in a trapped Bose-Einstein 
condensate with attractive interparticle interactions in 
the framework of mean-field theory (GP equation). We 
can conclude that, in general, only a fraction of the con- 
densate takes part in the collapse, while the rest of the 
particles are screened from the emerging singularity by 
an effective potential barrier due to 'quantum pressure' 
effects. This potential barrier is somewhat reminiscent 
of the 'event horizon' that appears around gravitational 
singularities. We could also show that, for reasonable 
initial widths of the condensate, the fraction of particles 
removed from the condensate as a result of collapse is 
an almost 'universal' number ~ 0.7. Even though we 
here derived this fraction of ejected particles only ap- 
proximately, our definition of the related particle number 
N e Eq. (|l4|), can be made mathematically rigorous. This 
number can be obtained from the solution of Eq. (|^) . The 
simplified model based on Gaussian trial wave function, 
which we proposed here, is not accurate if the 'coupling 
constant' go is in the vicinity of go cr - A more thorough 
analysis of Eq. (@) will be reported in a forthcoming pub- 
lication. 
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The final stage of the collapse takes place very fast, so a 
good numerical estimate for £* was already obtained from 
this seemingly rather crude approximation, which however 
allowed us to avoid the difficulties related to the mathematical 
singularity at t — i». 
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